Prompt and Afterglow Emission Properties of Gamma-Ray 
Bursts with Spectroscopically Identified Supernovae 
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ABSTRACT 



We present a detailed spectral analysis of the prompt and afterglow emis- 
sion of four nearby long-soft gamma-ray bursts (GRBs 980425, 030329, 031203, 
and 060218) that were spectroscopically found to be associated with type Ic su- 
pernovae, and compare them to the general GRB population. For each event, 
we investigate the spectral and luminosity evolution, and estimate the total en- 
ergy budget based upon broadband observations. The observational inventory 
for these events has become rich enough to allow estimates of their energy con- 
tent in relativistic and sub-relativistic form. The result is a global portrait of 
the effects of the physical processes responsible for producing long-soft GRBs. 
In particular, we find that the values of the energy released in mildly relativistic 
outflows appears to have a significantly smaller scatter than those found in highly 
relativistic ejecta. This is consistent with a picture in which the energy released 
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inside the progenitor star is roughly standard, while the fraction of that energy 
that ends up in highly relativistic ejecta outside the star can vary dramatically 
between different events. 

Subject headings: gamma rays: bursts - supernovae: general 

1. Introduction 

The discovery in 1998 of a Gamma-Ray Burst (GRB) in coincidence with a very unusual 
supernova (SN) of Type Ic (GRB 980425 - SN 1998bw; spectroscopically identified) was 
a turning point in the study of GRBs, offering compelling evidence that long-soft GRBs 
(Kouveliotou et al. 1993) are indeed associated with the deaths of massive stars (Galama 
et al. 1998). The large energy release inferred for the supernova suggested that GRBs are 
potentially associated with a novel class of explosions, having unusual properties in terms 
of their energy, asymmetry, and relativistic ejecta. More importantly, however, GRB 980425 
provided the first hint that GRBs might have an intrinsic broad range of energies: the total 
energy output in 7-rays (assuming an isotropic energy release) was only ~ 7 x 10 47 erg, 
some four orders of magnitude less energy than that associated with typical GRBs (Bloom, 
Frail & Kulkarni 2003). Finally, the fact that GRB 980425/SN 1998bw was located in a 
nearby galaxy with a redshift z = 0.0085 (Tinney et al. 1998; at 1 35.6 Mpc it remains the 
closest GRB to the Earth] gave rise to the possibility that such lower-energy bursts might be 
more common than had previously been thought, but harder to detect due to instrumental 
sensitivity. 

Unfortunately, during the elapsed eight years, very few SNe have been observed simul- 
taneously with GRBs. To date, three more nearby GRBs have been unambiguously, spectro- 
scopically identified with supernovae, two of which were discovered in 2003 (GRB 030329 - 
SN 2003dh; Stanek et al. 2003; Hjorth et al. 2003, GRB 031203 - SN 20031w; Tagliaferri et al. 
2003) and one in 2006 (GRB 060218 - SN2006aj; Pian et al. 2006). Each of these SNe is of the 
same unusual type as SN 1998bw. Note, however, that there is weaker photometric evidence 
that many other GRBs are accompanied by SNe, mainly by identification of a late time "SN 
bump" in the GRB optical afterglow lightcurve (e.g. Bloom et al. 1999; Galama et al. 2000; 
Zeh, Klose, & Hartmann 2004; Woosley & Bloom 2006). In this paper, we address only the 
four events with spectroscopically verified SN associations. Among these, only GRB 060218 
(interestingly at 143.2 Mpc for z = 0.0335, the closest after GRB 980425; Mirabal et al. 
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2006; Pian et al. 2006) had 7-ray energetics somewhat comparable to GRB 980425, and 
could possibly be added to the intrinsically-faint GRB sample. In the case of GRB 030329, 
the total energy release was at the low end of the typical range (-E 7i i S o ~ 10 52 erg), much 
higher than in the other three events. In fact, SN2003dh was obscured by the extreme op- 
tical brightness of the GRB afterglow and was only detected spectroscopically in the GRB 
optical lightcurve. Finally, the total energy release of GRB 031203 was intermediate between 
that of GRB 980425 and regular GRBs. Ramirez-Ruiz et al. (2005) argued that the faint 
GRB 031203 was a typical powerful GRB viewed at an angle slightly greater than about 
twice the half-opening angle of the central jet. 

In the present study, we consider the energy released during the GRBs, the afterglow, 
and the SN explosion for these four events in all wavebands, from 7-rays to radio waves, and 
we also estimate their kinetic energy content. The properties of the prompt and afterglow 
emission are described in detail in §§ 2 and 3, respectively, and are compared in § 4. The 
bolometric energy calculation and the evolution of the explosion responsible for their asso- 
ciated SNe are presented in § 5, while § 6 discusses the combined GRB-SN properties and 
their potential implications. We conclude in § 7 with a brief summary of our primary results 
and their implications. 

2. Prompt Emission 
2.1. GRB 980425 

GRB 980425 triggered the Large Area Detectors (LADs; 20-2000 keV) of the Burst 
and Transient Source Experiment (BATSE) on board the Compton Gamma-Ray Observa- 
tory, on 1998 April 25 at 21:49:09 UT (Kippen et al. 1998). The burst was also observed 
with the Gamma- Ray Burst Monitor (GRBM; 40-700 keV) and one of the two Wide Field 
Cameras (WFC; 2— 27keV) on board the BeppoSAX satellite. Figure 1 shows the bright- 
est BATSE/LAD and WFC lightcurves of the burst. Frontera et al. (2000) performed a 
broadband spectroscopic analysis of the event using the WFC and GRBM data. However, 
their analysis did not yield sufficiently constrained spectral parameters, due to the much 
lower sensitivity and poorer energy resolution of GRBM compared to those of the LADs. 
We performed here a broadband time-resolved spectral analysis of GRB 980425, combining 
the LAD and the WFC data, to better constrain the spectral parameters and their evolution 
during the burst. 

We used the High Energy Resolution Burst data of the brightest LAD for this event, 
with 128-channel energy resolution (Kaneko et al. 2006). This datatype starts after the 
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BATSE trigger time, which corresponds to the beginning of interval A, as shown in Figure 1. 
We modeled the LAD background with a low-order polynomial function on spectra taken 
~ 800 s before and ~ 200 s after the burst episode. The WFC-2 detected the burst at a 
fairly large off-axis angle of 18?7 with a light collecting area of only 10 cm 2 (6% of the on- 
axis value). However, since the WFC is an imaging telescope with a coded mask, all imaged 
data can be corrected for background without ambiguity. In addition to the background 
correction, the data were corrected for the detector deadtime as well as for the variations of 
the spectral response over the detector area. 

To assure good statistics, we binned both the LAD and WFC data into four time 
intervals (A, B, C, and D in Figure 1); two intervals each before and after the WFC lightcurve 
peak, similar to the ones used by Frontera et al. (2000). We analyzed the data from both 
detectors jointly, using the spectral analysis software RMFIT (Mallozzi, Preece & Briggs 
2006). Since no significant signal was found above 300 keV in all intervals (< 2.2cr detection in 
each interval), we obtained the best fits using the "Comptonization" photon model (Kaneko 
et al. 2006): 



Here, A is the amplitude in photons s _1 cm -2 keV -1 , a is the power-law index, and £ pea k is 
the photon energy (£) where vF v oc £ 2 f{£) peaks. We found no indication that interstellar 
absorption was required in the fit and the absorption was, therefore, not used. Finally, no 
normalization factor was needed between the datasets of the two detectors. 

The spectral fitting results are summarized in Table 1. Figure 2 displays the evolution 
°f £pcak and a; both parameters clearly change from hard to soft, in line with the common 
GRB trend (e.g., Ford et al. 1995; Crider et al. 1997). Moreover, the hardest part of the 
burst has an £ pea k of 175 ± 13 keV, which is within la of the peak of the £ pea k distribution for 
time-resolved BATSE GRB spectra (Kaneko et al. 2006). Further, we found the £ pca k of the 
duration-integrated spectrum to be 122 ± 17keV (la uncertainty), placing the event within 
2a from the peak of the distribution of BATSE GRBs (Kaneko et al. 2006). We find an 
isotropic-equivalent total energy emitted in 1 — 10000 keV oi E j iso = (9.29 ±0.35) x 10 47 erg, 
using a luminosity distance of d L = 35.6 Mpc (for z = 0.0085; Tinney et al. 1998). 

The total energy fluence values in the X-ray and 7-ray bands are Sx (2-30 keV) — 1-99 x 
10~ 6 erg cm -2 and (30-400 keV) = 3.40 x 10~ 6 erg cm -2 , respectively. The ratio of the two, 



log(S x /S 7 ) = -0.23, makes GRB 980425 an X-ray rich GRB, following the definition of 
Lamb, Donaghy, & Graziani (2004) and Sakamoto et al. (2005). 
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2.2. GRB 030329 

GRB 030329 was detected with the FREGATE detector (8-400 keV) on board the High- 
Energy Transient Explorer (HETE-2) on 2003 March 29 at 11:37:14 UT (Vanderspek et al. 
2004, also Figure 3). A detailed spectral analysis of the event was performed by Vanderspek 
et al. (2004), who reported an £ pea k = 70.2 ± 2.3 keV for the duration-integrated spectrum. 
Based on their spectral parameters and flux values, we estimate E J: - lso = 1.33 x 10 52 erg 
(1-10000 keV) using d L = 790.8 Mpc (z = 0.1685; Greiner et al. 2003; see also Table 4). 

2.3. GRB 031203 

GRB 031203 was detected with the IBIS/ISGRI detector (15-500 keV) on the INTE- 
GRAL satellite on 2003 December 3 at 22:01:28 UT (Gotz et al. 2003). We obtained the 
IBIS/ISGRI data from the INTEGRAL Science Data Center, and processed them with OSA 
5.0. Figure 4 displays the lightcurve of this event, for which we estimated a T 90i7 duration of 
37. Oil. 3 s. We extracted a background spectrum at ~ 300 seconds prior to the burst and a 
source spectrum encompassing the Tg 0)7 ± 10 s. In addition, we extracted four time-resolved 
spectra, with durations determined by requiring their signal-to-noise ratios to be above 35, 
for sufficient statistics. The four intervals are indicated in Figure 4. 

We used XSPEC vl2.2.1 for all spectral analyses of this event. The extracted integral 
source spectrum (17— 500 keV; effective exposure time of 45.9 s) is best described by a sin- 
gle power law of index of —1.71 ± 0.08, consistent with the results of Sazonov, Lutovinov 
& Sunyaev (2004). To constrain the £ pea k, Sazonov, Lutovinov & Sunyaev (2004) fitted 
the spectrum with the Band function (two power laws smoothly joined at a break energy; 
Band et al. 1993) with fixed high-energy photon index, and derived £ pca k > 190 keV (90% 
confidence level). Employing the same method, we found £ pea k > 71 keV (la) and 36keV 
(90%), and a low-energy photon index of a = —1.39 ± 0.41. These results are consistent 
with Sazonov, Lutovinov & Sunyaev (2004) since the longer integration time used here (45.9 s 
versus 22 s) includes the softer tail portion of the event, while the background contribution is 
practically negligible. Using the single power-law fit parameters we estimate Sx/S^ = 0.49 
(Table 4), which is in the middle of the range predicted by Sazonov, Lutovinov & Sunyaev 
(2004). Therefore, although the INTEGRAL data alone cannot confirm an X- Ray-Flash 
(XRF) nature for GRB 031203, as suggested by Watson et al. (2004), we find the event to 
be clearly X-ray rich. 

The isotropic-equivalent total emitted energy was estimated extrapolating the single 
power law fit between 1 — 10000 keV in the source rest frame (z = 0.105; Prochaska et al. 
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2004) leading to an £ 7jiso = (1.67lg;^) x 10 50 erg, using d L = 472.6 Mpc. However, since we 
assumed no spectral break up to 10 MeV, we consider this estimate to be an upper limit of 
-E^iso- All time-resolved spectra were also best fitted with a single power-law model; the fit 
results are shown in Table 2 and the power-law index evolution is presented in Figure 5. In 
accord with Sazonov, Lutovinov & Sunyaev (2004) and as illustrated in Figure 5, we found 
no significant spectral evolution within the burst, although spectral softening is suggested 
by the lightcurve of this event (Figure 4; see also §2.5.2). 



2.4. GRB 060218 

GRB 060218 was detected with the Burst Alert Telescope (BAT; 15-350 keV) on board 
the Swift satellite, on 2006 February 18, at 3:34:30 UT (Cusumano et al. 2006). The burst 
was exceptionally long, with a T 90)7 duration of 2100 ± 100 s (Campana et al. 2006). The 
available BAT data types were Event data up to ~ 300 s after the trigger time (T ) and 
Detector Plane Histogram (DPH) data, thereafter. The latter data type has a variable 
time resolution; when in survey mode, the nominal resolution is 300 s, but shorter exposures 
(with a minimum of 64 s) are available immediately after a burst trigger 2 . After the burst 
trigger, the satellite slewed immediately to the BAT burst location and the X-Ray Telescope 
(XRT; 0.2— lOkeV) started observing the burst around T + 160 s. The XRT data were 
in Windowed Timing (WT) mode from 161 to 2770 s after T . Figure 6 shows the BAT 
and XRT WT lightcurves of GRB 060218. Notice the much higher photon flux in the X-ray 
lightcurve, which clearly indicates the XRF nature of this event. In fact, the burst was 
detected with BAT by an image trigger, which typically displays no significant signals in the 
BAT lightcurve (Cusumano et al. 2006). 

We downloaded the publicly available BAT and XRT data and processed them with 
FTOOLS v6.0.4. For the present analysis, we used the BAT Event data for the first 150 
seconds of the burst, and the BAT DPH and XRT WT data for the rest of the burst. This 
combination of burst trigger and survey mode led to 14 time bins. The time bins are also 
shown in Figure 6. Before extracting spectra, we rebinned the DPH data with baterebin 
to the correct energy calibration edges. Then we created the mask-weighting map for the 
burst location using batmaskwtimg to extract the mask- weighted spectra from the rebinned 
DPH data. Finally, we extracted the spectra using batbinevt for both the Event and DPH 
data, and generated the corresponding detector response matrices with batdrmgen. For the 
XRT data, we obtained the screened WT event file and extracted spectra for the 13 different 



2 BAT Data Products v3.1 by H. Krimm (http://www.swift.ac.Uk/BAT_Data_Products_v3.l.pdf). 
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time intervals corresponding to the DPH time intervals. We used xselect and an extraction 
region of 50". For each spectrum, we generated the corresponding Auxiliary Response Files 
(ARF) with xrtmkarf using the default empirical ARF file as an input. The spectra were 
then grouped so that each energy channel contained a minimum of 20 counts, and the latest 
available response matrix (v007) was used for the analysis. No pile-up corrections were 
necessary for this event. 

We analyzed with XSPEC vl2.2.1 the BAT and XRT data jointly, except for the first 
time interval for which only BAT Event data were available. The BAT and XRT energy 
ranges used for the analysis were 15— 150 keV and 0.6— 9keV, respectively. In the majority 
of these intervals, an absorbed power law with an exponential high-energy cutoff was the 
best fit, albeit with a significant low energy excess. This excess was also identified by 
Campana et al. (2006) and was fitted with an additional blackbody (BB) component. We 
have, therefore, included the BB in all spectral fits where the XRT data were included. 

The BAT Event data (T — 8 to T + 140) spectrum is best fitted with a single power 
law with a high-energy cutoff, with spectral index of —0.87 ± 0.75 and £ pC ak of 24.9 ± 
6.0 keV (see also Table 3). The estimated energy fluence during this interval (0.5— 150 keV) 
is (6.38 ±4.00) x lO^ergcm- 2 . From T + 140 to T + 2734, we analyzed the BAT-XRT joint 
spectrum; the best fit was an absorbed Band function with an additional BB component of 
kT = 0.150±0.004 keV. The low and high energy spectral index and the peak photon energy 
values were a = -1.44±0.06, (3 = -2.54±0.07, and £ pcak = 4.67±1.15 keV, respectively. The 
absorption model used was wabs, with the best-fit value of Njj = (0.60 ± 0.02) x 10 22 cm -2 . 
The unabsorbed energy fluence (0.5— 150 keV) was (1.72 loire) x 10~ 5 ergcm -2 , where the 
errors were estimated from the absorbed flux. We derived an isotropic-equivalent total 
emitted energy of E^ tiso = (4.33 i?;^) x 10 49 erg in 1 — 10000 keV (source frame), using 
d L = 143.2 Mpc (z = 0.0335; Mirabal et al. 2006; Pian et al. 2006). 

The results of the time-resolved spectral fits are listed in detail in Table 3. We fitted an 
absorbed power law with a high-energy cutoff and a BB component in all BAT-XRT joint 
spectra but the last two, for which a single power law with BB provided adequate fits. The 
N H values were found to be constant (N H ~ 0.6 x 10 22 cm" 2 ) throughout the first 10 BAT- 
XRT joint time intervals, and were fixed to that value for the last three spectra, where their 
best-fit Nh values were higher resulting in excessive values of unabsorbed flux. Fixing Nh 
did not affect significantly the other parameters in the fits. Figure 7 displays the spectral 
parameter evolution of the 14 time-resolved spectra. We observed a hard-to-soft spectral 
evolution in the non-thermal spectra, while the BB temperature {kT) remained constant. 
On the average, the BB component contributes 0.13% of the source flux throughout the 
burst prompt emission duration. 
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2.5. Comparison 

The 7-ray properties of the four events are summarized in Table 4, and Figure 8 shows 
the evolution of their isotropic-equivalent 7-ray luminosity, L 7) i SO (in 2— 500 keV, source- 
frame energy). For GRB 060218, the contributions from the thermal component are shown 
separately from the L 7i i SO of the non-thermal component. We find that GRB 060218 had 
similar 7-ray luminosity as GRB 980425; both were a few orders of magnitude fainter than the 
other two events. However, since GRB 060218 lasted considerably longer than GRB 980425, 
its E^iso was ~ 50 times larger. 

We should emphasize here that we measure directly only the energy radiated in the 
direction of the Earth per second per steradian per logarithmic energy interval by a source 
at luminosity distance di- The apparent bolometric luminosity may be quite different from 
the true bolometric luminosity if the source is not isotropic. 

To compare the lightcurves of these events, we fitted a two-sided Gaussian function to 
the main pulses of the lightcurves. We found that the time profiles of these events are similar 
to the overall profile found for GRBs (Nemiroff et al. 1994), with the rising part having a 
Half Width Half Maximum (HWHM) much smaller (< 1/2) than the HWHM of the decaying 
part. 

2.5.1. Duration-Integrated Spectra 

Figure 9 compares the best-fit duration-integrated spectra (2 — 500 keV in source frame) 
of the four events in vF v , which shows the energy radiated per logarithmic photon energy 
interval. We plot here the unabsorbed spectral models in each case (solid lines) together with 
the spectrally deconvolved, blue-shifted data points shown in gray. The data are binned here 
for display purposes. No absorption correction was necessary in the combined LAD-WFC 
fit of GRB 980425. We used the published spectrum for GRB 030329, and we extrapolated 
our best-fit model for GRB 031203 to the lower energy range shown here. For GRB 060218 
the absorption correction was significant as seen from the comparison of the model to the 
XRT data below 10 keV. 

In Figure 10 we compare the vF v peak energy, £ pe ak, of the same duration- integrated 
spectra, to the £ pcak distributions of 251 bright BATSE GRBs (Kaneko et al. 2006) and 37 
HETE-2 GRB/XRFs (Sakamoto et al. 2005). The £ pcak values of the four events are lower 
(softer) than the average value of the bright BATSE GRBs but they all seem to fit well 
within the HETE-2 £ pea k distribution. Note that the BATSE £ pe ak distribution shown here 
is derived using the brightest BATSE GRBs, which tend to be spectrally harder than dim 
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GRBs (Mallozzi et al. 1995). 

Finally, in Figure 11, we show where these four events fall in the ^peak-^.iso plane, 
in comparison to the empirical correlation (dashed line) found by Amati et al. (2002) and 
Lloyd- Ronning & Ramirez- Ruiz (2002). GRB 980425 lies farthest from the so-called Amati 
relation. Although GRB 031203 also lies away from the relation, this deviation could be the 
result of a viewing angle slightly outside the edge of the jet (Ramirez- Ruiz et al. 2005). In 
fact, we find a smaller lower limit on £ pea k for this event than previously reported, so it is 
possible that, after correcting for an off-axis viewing angle, it would be consistent with the 
relation (compare our Figure 11 to the left panel in Figure 3 of Ramirez- Ruiz et al. 2005). We 
also note that for GRBs 980425 and 031203, Ghisellini et al. (2006) have recently suggested 
two other possible scenarios that may cause such a deviation from the Amati relation: one 
is that the GRB radiation may have passed through a scattering screen, and the other is 
that we may have missed softer emission from the GRBs lasting much longer. In both cases, 
observed spectra could have much higher £ pca k values than the intrinsic ones. Altogether, 
the four events do not appear to follow a specific pattern or correlation in the £ pea k--E 7 ,i SO 
plane. 



2.5.2. Spectral Evolution 

Figure 12 shows the hardness ratio evolution of all four events, in their source- frame time. 
The hardness ratio is defined here as the ratio of the energy flux between 50— 500 keV and 
2— 50keV in the source-frame energy. The hard-to-soft evolution is seen in all events except 
GRB 031203. The lack of spectral evolution in GRB 031203 is probably due to the fact that 
this burst was relatively dim and the observed energy range was fairly narrow. The burst was 
only detected below 200 keV (see Figure 4) and the break energy seemed to always lie around 
or above this energy; therefore, the flux above this energy may have been overestimated. We 
also show in Figure 13, the evolution of the £ pea k i n three events (GRBs 980425, 030329, 
and 060218) for which the £ pea k values were determined from their time-resolved spectra. To 
characterize the decaying behavior, we fitted a power law (£ pca k oc t^) to each of the data 
sets. Since we are only interested in the decaying behavior, the first points of GRBs 980425 
and 060218 were excluded from the fits. We found = -1.46 ± 0.12, -1.17 ± 0.08, and 
— 1.40 ± 0.06, respectively; the best-fit power laws are shown as dotted lines in Figure 13. 
We note that in the external shock model, for example, if £ pea k is identified with the typical 
synchrotron frequency u m , it is expected to scale as t~ 3 / 2 once most of the energy is transfered 
to the external medium and the self-similar deceleration phase sets in (Blandford & McKee 
1976). This decaying index is actually similar to what we observed here, especially when 
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taking into account that the asymptotic power law of t 3 / 2 is expected to be approached 
gradually rather than immediately. 

3. Afterglow Emission 

3.1. X-ray Afterglow 

Figure 14 shows the evolution of the prompt and afterglow isotropic-equivalent lumi- 
nosity in 0.3— lOkeV (Xx,iso) i n the source-frame time for all four events. We estimated the 
prompt X-ray luminosity values from the spectral analysis of the prompt 7-ray emission (§2) 
by extrapolating the best-fit 7-ray spectra to the X-ray energy range for each event. The 
afterglow X-ray luminosity values were obtained in various ways. 

For GRB 030329, we derived the values using the data presented in Tiengo et al. (2003, 
2004). In the case of GRB 031203, we analyzed the (now) archival XMM and Chandra 
observations. Our analysis results are consistent with those presented in Watson et al. 
(2004). For GRB 060218, we analyzed the Swift XRT photon counting (PC) mode data, 
starting ~ 6200 s after the burst trigger. We extracted two time-averaged spectra, covering 
6191— 8529 s and 1.1 x 10 4 — 1.3 x 10 5 s after the burst trigger, respectively, to account 
for spectral evolution, and fitted each spectrum with a single power law. The spectral 
parameters were then used to convert the count rates to the luminosity values. The PC 
data after 1.3 x 10 5 s were associated with relatively large uncertainties due to the fact that 
the count rates were approaching the XRT sensitivity limit, and therefore, we used the 
fitted parameters of the second spectrum above to estimate the luminosity values for these 
data points. We also include in this plot the two Chandra observations of GRB 060218 at 
~ 10 6 s after the burst trigger, which were presented in Soderberg et al. (2006). Finally, 
we reanalyzed the BeppoSAX observations of GRB980425/SN 1998bw, which were initially 
presented in Pian et al. (2000). We used only the data from the BeppoS AX/Medium Energy 
Concentrator Spectrometers and followed closely the procedure described in Pian et al. 
(2000). The only potentially significant differences between the current and the analysis 
of Pian et al. (2000) is that we estimated the background using the same observation, while 
Pian et al. (2000) used calibration blank-sky observations. As a result, we found slightly 
lower net count rates (1.6— 10.0 keV) during the two observations in April 1998 (as compared 
to Table 1 in Pian et al. 2000), indicating a slightly flatter decay than was initially reported. 
In addition, we found that the observation beginning April 27 1998 contained 6.1 ks less data 
than previously reported (total of 14,862 s exposure time). The count rates and exposure 
times for the remaining observations are consistent with Pian et al. (2000). 
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To investigate the energy dissipation behavior in the X-ray afterglow, we fitted a natural 
cubic spline function to the Lx,iso history for each individual afterglow and estimated the 
cumulative emitted energy as a function of time. The evolution of the cumulative X-ray 
afterglow energy for all four events is shown in Figure 15. The integration time intervals 
varied from event to event, since the data time coverage was different for each event. We 
defined the integration time (i.e., "afterglow" time interval) of each event to be from the 
earliest X-ray observation time (~ 10 4 s) to the time of the latest available data point. In the 
case of GRB 060218, however, it was difficult to determine the starting point of the afterglow 
since the prompt emission was extremely soft and long and the X-ray follow-up started much 
earlier than in the other three events. It is possible that the steep decay observed before 
~ 10 4 s (see Figure 14) belongs to the decaying part of the prompt emission; we, therefore, 
chose the afterglow starting time for this event to be 1.1 x 10 4 s, which corresponds to the 
onset of the break in its Lx,iso history. 

To compare the cumulative energy evolution to regular GRBs, we also show in Fig- 
ure 15 the cumulative energy evolution of the 10 Swift GRBs with known redshifts that 
were presented in Nousek et al. (2006). For these 10 events, we converted the Lx,i SO values 
in 2— lOkeV presented there, to the Lx.iso m the energy range we used for our four events 
here (i.e., 0.3— lOkeV), using the same spectral parameters used in Nousek et al. (2006), 
and again, fitted a spline function to each Lx.iso evolution. The start and end times of the 
integration were the first and the last points of the actual observations. Figure 15 shows 
that the final ^x.iso values for most of the 10 Nousek et al. (2006) GRBs span the range of 
~ 10 50 - 10 52 erg. Out of the four SN-GRB events, only GRB 030329 falls within this range, 
while the other three events fall between ~ 10 48 and ~ 10 49 erg. We note that there exists 
a selection effect based on the observed photon flux: an event would be more likely to be 
detected when it is closer to us than farther, for a given intrinsic luminosity. Therefore, the 
difference could be partially due to the fact that these four events occured relatively nearby 
compared to the 10 Swift events. From the cumulative X-ray energy evolution, we also de- 
rived T 9 o,x values for all events, which we define here as the time in the source frame during 
which 90% of the radiated afterglow energy (0.3— lOkeV in the source frame) is accumulated. 
Table 4 displays the comparison of the X-ray afterglow properties of all four events, along 
with their 7-ray prompt properties. Again for the purpose of comparison, Figure 16 shows 
the -E 7i iso, Ex,\so, and Tgo,x values of the four SN-GRBs, along with those of the 10 Nousek 
et al. (2006) GRBs. Here we plot the -E 7i i SO values between 20— 2000 keV (in the source 
frame). There is a clear correlation between E^ : i SO and -Ex.iso (E lt i so <x E^ ^' 3 ), and an 
anti-correlation between -E 7) i SO and T 90i x (-E^iso oc ^go'x^ 3 )- We also find T 90j x oc E^'J^ ' 2 , 
which can be interpreted as follows. 

Events that have large isotropic-equivalent energy (both in 7-rays, E lt i so , and in the 
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kinetic energy of their afterglow, -Ek,iso) have a large Ex,i SO , indicating a reasonably narrow 
spread in the efficiency of converting the afterglow kinetic energy into radiation. Moreover, 
they are typically associated with narrow jets. This means that most of their kinetic energy 
is in relativistic outflow carried by highly-relativistic ejecta (with r ~ 30 — 50; Granot & 
Kumar 2006). Therefore, they have a relatively small Tgo,x; the time at which most of 
the energy is injected into the afterglow shock, of the order of 10 3 s < tbrcak,2 < 10 4 s 
(Nousek et al. 2006). On the other hand, events that have a small -Ek,iso and small E J:iso , 
naturally have a small £x,i so - Such events also tend to have most of their relativistic outflow 
energy residing in mildly relativistic ejecta, rather than in highly relativistic ejecta (Granot 
& Ramirez- Ruiz 2004; Waxman 2004). As a result, most of the kinetic energy is transfered 
to the afterglow shock at relatively late times, hence a large Tgo,x- Further, since such events 
tend to be only mildly collimated, the true total energy in relativistic (r > 2) ejecta has a 
significantly smaller spread than the one of the isotropic equivalent energy (E J>iso or Ex,i so )- 



3.2. Radio Afterglow 

We performed radio observations of GRB 060218 with the Westerbork Synthesis Ra- 
dio Telescope (WSRT); these results are tabulated in Table 5. Additional radio data for 
this burst and radio data for the three other events were collected from the literature 
(GRB 980425; Frail et al. 2003 and Kulkarni et al. 1998, GRB 030329; Berger et al. 2003, Frail 
et al. 2005, Taylor et al. 2005, van der Horst et al. 2005, and Resmi et al. 2005, GRB 031203; 
Soderberg et al. 2004, GRB 060218; Soderberg et al. 2006 and Kamble et al. 2006). The radio 
flux lightcurves in 4.8/4.9 GHz, at which all of the four events were observed, are shown in 
Figure 17. 

As with the X-ray afterglow, we fitted the flux lightcurve of each event (in a wide band, 
where available) with a natural cubic spline, to estimate the cumulative energy evolution 
emitted in radio. To minimize the effect of scintillation in the radio lightcurves, we used 
fewer nodes than data points in the fitting. This resulted in a smooth fit to the data that 
retains the overall evolution of the lightcurve. The curve is forced to be at very early 
times (< 0.01 days) and late times (> 500 days), or to the inferred flux of the host galaxy 
in the case of GRB 031203 (Soderberg et al. 2004). The fitted lightcurves were then used to 
create a grid in time-frequency space to obtain the flux profile in a wider frequency range 
(Figure 18). Using this, we constructed the cumulative radio afterglow energy evolution 
in 5— 7 GHz, shown in Figure 19. In Table 4, we also present the comparison of the radio 
afterglow properties of the four events, along with their 7-ray prompt and X-ray afterglow 
properties. Similar to T 9 o,x defined above, we define T 9 o,r here as the time in the source 
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frame during which 90% of the radiated afterglow energy is accumulated in 5—7 GHz in the 
source-frame energy. The integration time of the radio afterglow emission is 1 day to 100 days 
after the burst trigger time, in the rest frame of the source. 

As can be seen in Table 4, the isotropic equivalent energy that is radiated in the radio 
(£ , R ) i SO ) is ~ 3.5 — 5 orders of magnitude smaller than that in X-rays, Ex,i SO - This is pre- 
dominantly due to the fact that vF v typically peaks closer to the X-rays than to the radio, 
and it is very flat above its peak while it falls much faster toward lower energies. Another 
effect that enhances the difference between -E^iso and Ex,i SO is that the latter is calculated 
over a much wider energy range (0.3 — 10 keV versus 5 — 7 GHz). Finally, since these are 
isotropic equivalent energies, most of the contribution to -Er^o is from significantly later 
times than for Ex,i SO , and the collimation of the outflow together with relativistic beam- 
ing effects could result in much larger -Ex,iso than E Rjiso . We note that for at least two 
(GRB 908425 and GRB 060218) of the SN-GRBs, the isotropic-equivalent emitted energy at 
optical wavelengths can be (much) larger than Ex,iso because it is typically dominated by the 
contribution from the SNe (powered by radioactive decay) rather than by the GRB afterglow 
emission (see also § 5). 



4. Prompt and Afterglow Properties 

The isotropic-equivalent luminosity of GRB X-ray afterglows scaled to t = 10 hr after 
the burst in the source frame, Lx,i SO (10 hr), can be used as an approximate estimator for 
the energy in the afterglow shock for the following reasons. First, at 10 hr the X-ray band 
is typically above the two characteristic synchrotron frequencies, v m (of the accelerated 
electron with minimum energy) and v c (of the electron whose radiative cooling time equals 
the dynamical time), so that the flux has very weak dependence 3 on e# and no dependence on 
the external density, both of which are associated with relatively large uncertainties. Second, 
at 10 hr the Lorentz factor of the afterglow shock is sufficiently small (r ~ 10) so that a 
large fraction of the jet is visible (out to an angle of ~ ~ 0.1 rad around the line of 
sight) and local inhomogeneities on small angular scales are averaged out. Finally, the fact 
that the ratio of Lx,i SO (10hr) and E^^ Q is fairly constant for most GRBs, suggests that both 
can serve as a reasonable measure of the isotropic-equivalent energy content of the ejected 
outflow. 

Figure 20 shows Lx,iso(10hr) in 2—10 keV rest-frame energy as a function of their 



3 When synchrotron-self Compton is taken into account the dependence on eg becomes much stronger 
(Granot, Konigl & Piran 2006). 
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isotropic 7-ray energy release (-E 7 ,i SO , in 20—2000 keV) for the four events, in comparison to 
regular GRBs. A linear relation, L x ,iso(10 hr) oc E^^ iso , seems to be broadly consistent with 
the data, probably suggesting a roughly universal efficiency for converting kinetic energy 
into gamma rays in the prompt emission for these four events and those of "regular" long 
GRBs. This "universal" efficiency is also likely to be high (i.e., the remaining kinetic energy 
is comparable to, or even smaller than, the energy dissipated and radiated in the prompt 
emission). If this is the case, the well-known efficiency problem for long GRBs also persists 
for SN-GRB events. Not surprisingly, there is also a correlation between the cumulative 
X-ray afterglow energy and isotropic 7-ray energy (see Figure 16 and also § 6). 



5. SN Properties 

Pian et al. (2006) presented the bolometric luminosity evolution (in optical and near- 
infrared wavelengths; 3000— 24000 A) for all the SNe associated with the GRBs considered 
here, namely SNe 1998bw, 2003dh, 20031w, and 2006aj. We fitted a cubic spline to the 
luminosity evolution of each SN (Figure 3 of Pian et al. 2006), and estimated the cumulative 
energy emitted by the SN as a function of source-frame time. To prevent the spline fits to 
diverge after the last real data points (< 40 days), we assumed that the luminosity becomes 
much smaller (<C 10 41 erg s^ 1 ) at much later time (^> 1000 days). The cumulative SN 
energy history is plotted in Figure 21 for all four SNe. We also estimated Tg^sN, during 
which 90% of the SN energy is accumulated in 3000— 24000 A since 1 day after the burst, 
as well as the total energy emitted in the same energy range, -EsN.iso- These values are 
presented in Table 4, along with the burst properties in other wavelengths. It should be 
noted here that in most cases, -EsN.iso is monopolized by the SN emission while the GRB 
afterglow dominates at all other wavelengths. The integration time for the SN cumulative 
energy was from 1 day to 100 days after the burst trigger. The cumulative energy of the four 
SNe span a much narrower range than the one reported for the entire SNe Ic class (Mazzali 
et al. 2006a). Among these four, SN2006aj is the least energetic, albeit more energetic than 
other broad-lined SNe and normal SNe (Pian et al. 2006). 

Furthermore, from the radio observations of type Ib/c SNe presented by Weiler et al. 
(2002), we note that the peak radio luminosity of SN1998bw observed at 6cm (5 GHz) was 
a few orders of magnitude higher than the other type Ib/c SNe, although the overall time 
evolution of this event was comparable to the others. 
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6. Discussion 

The compilation of the radiated energy inventory, presented in Table 4 and Figure 22, 
offers an overview of the integrated effects of the energy transfers involved in all the physical 
processes of long GRB evolution, operating on scales ranging from AU to parsec lengths. 
The compilation also offers a way to assess how well we understand the physics of GRBs, by 
the degree of consistency among related entries. We present here our choices for the energy 
transfers responsible for the various entries in the inventory. Many of the arguments in this 
exercise are updated versions of what is in the literature, but the present contribution is the 
considerable range of consistency checks, which demonstrates that many of the entries in the 
inventory are meaningful and believable. 

In the following, § 6.1 repeats the minimal Lorentz factor estimates in Lithwick & Sari 
(2001) using the parameters inferred for all of the four SN-GRB events in this work. § 6.2 
compares a simple model for the emitted radiation to the observational constraints on the 
integrated energy radiated in the afterglow phase. Finally, the energy contents in relativistic 
and sub-relativistic form are estimated in §§ 6.3 and 6.4, respectively. 

6.1. Constraints on the Lorentz Factor 

As is well known, the requirement that GRBs are optically thin to high-energy photons 
yields a lower limit on the Lorentz factor of the expansion. We estimate here the minimal 
Lorentz factor of the outflow based on our observational analysis, by using the results of 
Lithwick & Sari (2001). They derived two limits; A and B. Limit A is from the requirement 
that the optical depth to pair production in the source for the photon with the highest 
observed energy, £ max , should be smaller than unity. Limit B is from the requirement that 
the Thompson optical depth of the e ± pairs that are produced by the high-energy photons 
in the source does not exceed unity. The required parameters to estimate these limits are 
the variability time of the 7-ray lightcurve, ST (which is related to the radius of emission 
by R ~ T 2 cST), and the photon flux per unit energy, dN p ^/ dAdtdSp^ = fS'^ for £ min < 
£ p h < £ m ax : i n the latter, both the normalization / and the high energy photon index f3 
are needed. Although the result depends on the exact choice of parameters, representative 
values are presented in Table 6. For each event, we determined ST to be the Full Width Half 
Maximum (FWHM) of a two-sided Gaussian function fitted to the prompt 7-ray lightcurve, 
the peak of which aligns with the main peak of the lightcurve. We then analyzed the "peak" 
spectrum with the integration time of ST to obtain the high-energy photon index (5 and the 
£ P eak values shown here. The minimal values of the Lorentz factor derived in this way are 
generally smaller for less energetic events (in terms of E^ :iso or £x,iso)- F° r the event with 



-16- 



the smallest E~ /jiso , GRB 980425, we find r min ~ 2.4 if we adopt our lower limit for the high- 
energy photon index ((3 > 3.5; la), while a larger value of j3 would lower r min (consistent 
with the results of Lithwick & Sari 2001). 



6.2. Integrated Radiated Energy 

The integrated (isotropic-equivalent) radiated energy during the afterglow, is given by 
-Erad,iso = J dt J dv L v - mo {t) 

^t^peak-^iso (^peak) ^t^i/^peak^peak(^peak)-^- / ;/,iso(^peak) ■> (1) 

where t, u, and L uiso are measured in the cosmological frame of the GRB, t pC ak is the time 
when tL iso (t) peaks, v is the frequency where f-^iso^peak) peaks; 

^ = [dt L iso (t) ^ 

^peak-^-'iso(^peak) 

and 

^ -^iso(^peak) 

^peak (^pcak ) -^ivso (^pcak ) 

are factors of order unity here (although we typically expect r\ t r\ v ~ 10). Since in practice 
vL^sso peaks near the X-ray band, we can assume that Lx,iso only mildly underestimates 4 L iso , 
and the time t pea k when tL; so (t) peaks is usually rather close to the time tx when tLx,i SO (t) 
peaks. Therefore, ^x-^x,iso(^x) provides a convenient lower limit for -E ra d,iso within an order 
of magnitude, although it is still possible that the X-ray observations might have missed the 
actual time when most of the energy was radiated, resulting in a significant underestimate. 
The values of tx-^x,iso(^x) are provided in Table 7. We find that these values are typically a 
factor of ~ 2 — 3 smaller than Ex,i SO , suggesting r\ t ~ 2 — 3. Although rj t is defined for the 
bolometric luminosity rather than for the X-ray luminosity, the values derived here are still 
fairly representative. 



6.3. Energy Inventory - Relativistic Form 

Supernova remnants are understood reasonably well, despite continuing uncertainty 
about the initiating explosion; likewise, we hope to understand the afterglow of GRBs, 



4 Evcn if vF v peaks below the X-rays, it is very flat above its peak, so a significant fraction of the afterglow 
energy is still radiated in the X-rays. 
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despite the uncertainties about their trigger mechanism. The simplest hypothesis is that the 
afterglow is due to a relativistic expanding blast wave. 5 The complex time structure of some 
bursts suggests that the central trigger may continue (i.e., the central engine may remain 
active) for up to ~ 100 seconds. However, at much later times all memory of the initial 
time structure would be lost. All that matters then is essentially how much energy has been 
injected and its distribution in angle and in Lorentz factor, e(> T, 9), where e = dE/dQ. 

6.3.1. Kinetic Energy Content 

An accurate estimate of the kinetic energy in the afterglow shock requires detailed 
afterglow modeling and good broadband monitoring, which enables one to determine the 
values of the shock microphysical parameters (electron and magnetic energy equipartition 
fractions, e e and e#, and shock- accelerated electron power-law index, p). However, even 
then, it provides only a lower limit for the kinetic energy due to the conventional and highly 
uncertain assumption that all of the electrons are accelerated to relativistic energies (Eichler 
& Waxman 2005; Granot, Konigl & Piran 2006). Nevertheless, an approximate lower limit on 
the isotropic-equivalent kinetic energy in the afterglow shock, -Ek,iso, can be obtained from 
the isotropic-equivalent X-ray luminosity, Lx t i SO , since the typical efficiency of the X-ray 
afterglow, ex = tLx t i SO {t) / E kiso (t) , is < 10~ 2 (Granot, Konigl & Piran 2006): a rough lower 
limit on E^ [so is obtained by adopting e(tx) ~ 10~ 2 , where tx is defined in the previous section 
and the values for the four events are presented in Table 7. For GRBs 980425 and 060218, 
we estimate E k: - lso (tx) ~ 6 x 10 49 ex j 1 _ 2 er g> an d E k: - lso (tx) ~ 3 x 10 49 ex j 1 _ 2 erg, respectively, 
where ex -i = ex(tx)/10~ 2 . These rough estimates are similar to the energy inferred from a 
more detailed analysis of the X-ray and radio observations of this event (-Ek,iso ~ 5 x 10 49 erg; 
Waxman 2004). For GRB 031203, we find E kiiso (t x ) ~ 2 x 10 50 e x ^_ 2 erg. 

For GRB 030329, we derive E k>iso (tx) ~ 5 x 10 52 e x 1 _ 2 erg. This estimate is comparable 
to that from the broadband spectrum at t !=s 10 days (E k: - lso ~ 5 x 10 52 erg for a uniform 
density and ~ 10 52 erg for a wind; Granot, Ramirez- Ruiz & Loeb 2005), assuming negligible 
lateral expansion of the jet (Gorosabel et al. 2006). For rapid lateral expansion, the inferred 
value of i?k,iso(10 days) is lower (~ 1.6 x 10 51 — 5 x 10 51 erg, Berger et al. 2003; Granot, 
Ramirez- Ruiz & Loeb 2005) but should correspond to a similar -Ek.iso before the jet break time 
(tj 0.5 days), for a comparable initial half-opening angle. There were, however, several 



5 For GRB 030329 this picture is supported by direct measurements of the angular size of its radio afterglow 
image (Taylor et al. 2004, 2005), which show a superluminal apparent expansion velocity that decreases with 
time, in good agreement with the predictions of afterglow models (Oren, Nakar & Piran 2004; Granot, 
Ramirez- Ruiz & Loeb 2005). 
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re-brightening episodes observed in the optical afterglow lightcurve of GRB 030329, between 
~ 1.5 days and a week after the GRB (Lipkin et al. 2004), suggesting energy injection into 
the afterglow shock that increased its energy by a factor of ~ 10 (Granot, Nakar & Piran 
2003). This would imply that -Ek,iso at tx ~ 2 x 10 4 s was a factor of ~ 10 lower than our 
rough estimate, or that €x(tx) is as high as ~ 0.1. A possible alternative explanation for 
the relatively high value of tx-^x,iso(^x) which does not require a high afterglow efficiency 
(ex) comes about if a flare in the X-rays was present around tx, which was not detected 
in the optical lightcurve at the similar time. In this case, I/x,iso(^x) would be dominated 
by late time activity of the central source rather than by emission from the external shock 
(Ramirez- Ruiz et al. 2001; Granot, Nakar & Piran 2003). 



6.3.2. Minimal Energy Estimates 

Here we derive a simple but rather robust estimate for the minimal combined energy 
in the magnetic field (E B ) and in the relativistic electrons (E e ) that are responsible for the 
observed synchrotron emission of flux density F Vj r at some frequency ur, E m i n = min(i?B + 
E e ). This estimate is applied for the sub-relativistic flow, to avoid the effects of relativistic 
beaming and reduce the uncertainty on the geometry of the emitting region. Since the source 
is not resolved, the energy is estimated near the time of the non-relativistic transition, t-^n, 
where we have a handle on the source size. 

We follow standard equipartition arguments (Pacholczyk 1970; Scott & Readhead 1977; 
Gaensler et al. 2005; Nakar, Piran & Sari 2005). The minimal energy is obtained close to 
equipartition, when Es/E e = 3/4. At such late times (t ~ £nr ^ 10 2 days) it is easiest 
to detect the afterglow emission in the radio, so vr would typically be in the radio band. 
Furthermore, ^(^nr) is also usually around the radio band, meaning that the electrons 
radiating in the radio would carry a reasonable fraction of the total energy in relativistic 
electrons. Still, the total energy of all electrons would be a factor of > 10 larger than E e , 
which would increase the total minimal energy by a factor of > 6. Since the kinetic energy is 
expected to be at least comparable to that in the relativistic electrons and in the magnetic 
field, the total energy is likely to be at least an order of magnitude larger than E m - m . 

Following Nakar, Piran & Sari (2005), E e = N e ^ e m e c 1 ^ where N e is the number of 
electrons with a synchrotron frequency z/ syn ~ v R , and therefore, with a Lorentz factor 
7 e Ri [2irm e c(l + z)u R /eB] 1/2 . Also, since F VjR ps N e P U)inax (l + z)/4ird 2 L where P^ max ~ 
-Psyn(7e)/t / syn(7e) ~ o- T m e c 2 B /3e, N e Ylnd 2 L eF v> R/ {1 + z)<7Tm e c 2 B. Estimating the emit- 
ting volume as V = (4n/3)R 3 /i] with 77 = IO771 (rj' 1 is the fraction of the volume of a sphere 
with the radius R of the emitting material that is actually occupied by the emitting mate- 
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rial), while in terms of our observed quantity, £nr, -R(^nr) = oc^nr/ (1 + z) with a = 2ao.3 (a 
is the average apparent expansion velocity at £nr m the cosmological frame of the GRB/SN 
in units of the speed of light), we obtain 

The resulting estimates of E m[n for the four events discussed in this paper are given in 
Table 8, along with £nr, and F v ^ values. We used ur = 4.86 GHz. As a sanity check, 
we calculate the minimal external density that corresponds to a total energy of lO-E^n; 
^min = 10-E m in/(47r/3)i?NR m p c2 > where i?NR = -R(^nr) = octNR/(l + z )-> m which we use 
our fiducial values of a = 2 and rj = 10. These values are also presented in Table 8. It 
is useful to compare E m - m to other energy estimates. For GRB 980425, the X-ray afterglow 
observations suggest an energy of ~ 5 x 10 49 erg ~ 30-E min in a mildly relativistic roughly 
spherical component (Waxman 2004). For GRB 030329 the total kinetic energy at late 
times is estimated to be ~ 3.2 x 10 50 erg ~ 100i? m i n (Granot, Ramirez- Ruiz & Loeb 2005). 
Therefore, it can be seen that while the ratio of the total kinetic energy in relativistic outflow 
(7/3 > 1) around £nr and E min is indeed > 10, it is typically < 100. Moreover, the fact that 
these different energy estimates are consistent lends some credence to these models. 

Some cautionary remarks are in order. The above calculation is only sketchy and should 
be taken as an order of magnitude estimate at present. For example, the usual assumption 
that at £nr the flow is already reasonably well described by the Newtonian spherical Sedov- 
Taylor self-similar solution, for which r\ fa 10, is probably not a very good approximation. 
Numerical studies show that there is very little lateral expansion of the GRB jets while 
the flow is relativistic, and therefore it takes at least several dynamical timescales after t NR 
for the flow to approach spherical symmetry. Furthermore, since the flow is still mildly 
relativistic at this stage, there is still non-negligible relativistic beaming of the radiation 
toward the observer from the forward jet and away from the observer from the counter-jet. 
Altogether, the flux is somewhat enhanced due to this mild beaming, and the fraction of the 
total solid angle that is occupied by the flow at t NR is still considerably smaller than unity. 
Consequently, this introduces an uncertainty of at least a factor of a few in the estimate of 
E m in', however, since the total energy is expected to be > 10E mm , it would still be at least a 
factor of a few larger than our estimate of E m i n . The theoretical uncertainty on the dynamics 
of the flow should improve with time as more detailed numerical studies become available. 

Another important uncertainty is in the determination of the non-relativistic transition 
time, £nr- The better estimate of E min should become obtainable as more well-sampled 
afterglow observations are made, and the modeling gets more precise so that one can more 
carefully estimate both t NR and F^ R (t NR ). For GRB 980425 we related t NR to the time at 
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which the X-ray lightcurve steepens, which likely corresponds to the deceleration time of 
the mildly relativistic ejecta. For GRB 030329, £nr was selected to roughly correspond to 
the time at which the radio lightcurve flattens and at the same time to be consistent with 
estimates derived from direct size measurements of the event (Oren, Nakar & Piran 2004; 
Granot, Ramirez- Ruiz & Loeb 2005). In the other two cases, only a crude estimate of £nr 
can be made as it is unclear whether current observations clearly display a signature of the 
non-relativistic transition. 



6.4. Energy Inventory - Sub-relativistic Form 

Despite the wide range in energies in relativistic ejecta, and even wider range in E ltiso , 
the total (non-neutrino) energy of the associated SN in all four events spans, at most, a factor 
of 10. Most of this energy is in non-relativistic (7/5 < 2) kinetic energy: the integrated light 
of the SN is negligible. While not standard candles, the optical luminosities of the four SNe 
at peak are all much brighter than average Type lb or Ic SNe (Woosley & Bloom 2006; Pian 
et al. 2006; Ferrero et al. 2006). Since the brightness of Type I supernovae at peak is given 
by the instantaneous rate of decay of 56 Ni, the 56 Ni masses are thus inferred to be in the 
range of ~ (0.2 - O.7)M . 

To produce this much 56 Ni in a Wolf-Rayet (WR) star requires a kinetic energy of at 
least ~2x 10 51 erg, even in lower mass WR stars (Ensman & Woosley 1988). In higher mass 
stars, a still greater energy is required for the lightcurve to peak within two weeks after the 
maximum. Large kinetic energies are also inferred from detailed models of the explosion, 
especially the lightcurve and the velocity histories of spectral features (see Woosley & Bloom 
2006, and references therein). In summary, the supernova kinetic energies in the four well- 
studied events almost certainly lie within the range 2 x 10 51 — 2 x 10 52 erg. In fact, the range 
of typical GRB-SNe may be much smaller, with the brightness at peak varying by no more 
than one magnitude in all four events and the kinetic energy in at least three of the four 
events (all but SN2006aj) within a factor of two of 1.5 x 10 52 erg (Woosley & Bloom 2006). 

Finally, in Figure 23, we compare the collimation-corrected total emitted energy in 7- 
ray (E 7 ) and supernova kinetic energy estimates (£k) of these four SN-GRB events with 
regular GRBs and other broad-lined (1998-bw like) SNe without GRBs. For GRB 030329, 
we used the jet angle estimate of 9j ~ 0.083 — 0.14 rad (Gorosabel et al. 2006). For the other 
three events, the isotropic-equivalent emitted energy, E 7t i so , was used as an upper limit, as 
there was no observational evidence of jet breaks for these events. The E 1 values for 27 
regular GRBs were adopted from Ghirlanda, Ghisellini, & Lazzati (2004), where the -E 7i i S o 
was again used as an upper limit for GRBs with no 9j constraints. All the E k estimates 
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were taken from the literature (SN1998bw; Woosley, Eastman, & Schmidt 1999, SN2003dh; 
Mazzali et al. 2003, SN20031w; Deng et al. 2005, SN2006aj; Mazzali et al. 2006b, SN2002ap; 
Mazzali et al. 2002, SN2003jd; Mazzali et al. 2005, SN1997dq & SN 1997dq; Mazzali et al. 
2004, and SN2005bf; Tominaga et al. 2005; Folatelli et al. 2006). We note that the explosion 
energy (i.e., SN E^) is much larger than the energy released as GRBs (E 7 ), and spans much 
narrower range than E 1 . 



7. Concluding Remarks 

One of the liveliest debated issues associated with GRBs is on the total energy released 
during the burster explosion: are GRBs standard candles? The GRB community has vacil- 
lated between initial claims that the GRB intrinsic luminosity distribution was very narrow 
(Horack et al. 1994), to discounting all standard candle claims, to accepting a standard total 
GRB energy of ~ 10 51 ergs (Frail et al. 2001), and to diversifying GRBs into "normal" and 
"sub-energetic" classes. The important new development is that we now have significant 
observational support for the existence of a sub-energetic population based on the different 
amounts of relativistic energy released during the initial explosion. A network of theoretical 
tests lends credence to this idea. The existence of a wide range of intrinsic energies that we 
presented in this work may pose challenges to using GRBs as standard candles - it is also 
worth stating explicitly that, when viewed together, these four events fall away from the 
Amati relation. 

Our results are consistent with the emerging hypothesis that GRBs and XRFs share a 
common origin in massive WR stars. The central engine gives rise to a polar outflow with 
two components (Woosley & Bloom 2006). One large angle outflow (the SN), containing 
most of the energy and mass, is responsible for exploding the star and producing the 56 Ni 
to make the SN bright. Only a tiny fraction of the material in this component reaches 
mildly relativistic velocities, which is more narrowly focused. A second outflow component 
(the GRB jet) occupies a narrower solid angle, probably contains smaller energy (which can 
range from comparable to much smaller), and most of its energy is in material with relativistic 
velocities (where the typical Lorentz factor of the material that carries most of the energy in 
this component can vary significantly between different SN-GRBs). After it exits the star, 
internal shocks within this jet and external shocks with the residual wind material around 
the star make the GRB or XRF and its afterglow. Apparently, the properties of the broad 
component are not nearly so diverse as those of the core jet (Ramirez-Ruiz & Madau 2004; 
Soderberg et al. 2006; Woosley & Bloom 2006). 

We have argued, using well-known arguments connected with parameters such as opacity 
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and variability timescales, that these less-energetic events do not require a highly-relativistic 
outflow. Our best estimates of Lorentz factors, T, for these events are in the range of 2—10. 
Indeed, it is much more difficult to produce a jet with very high Lorentz factor - i.e., a 
high energy loading per baryon - than with low Lorentz factor. A jet with low Lorentz 
factor could result even if a jet of relatively pure energy is produced, since it may be loaded 
with excess baryons by instabilities at its walls as it passes through the star, or if it does 
not precisely maintain its orientation (Ramirez- Ruiz, Celotti & Rees 2002; Aloy et al. 2002; 
Zhang, Woosley & Heger 2004). The above suggest that GRBs made by jets with lower 
Lorentz factor should be quite common in the universe (Granot & Ramirez- Ruiz 2004). 

Continued advances in the observations will surely yield unexpected revisions and addi- 
tions in our understanding of GRBs in connection with SNe: currently, we are attempting to 
draw large conclusions from limited observations of exceedingly complex phenomena. How- 
ever, the big surprise at the moment is that these SN-GRB events appear to be intrinsically 
different from and much more frequent (Granot & Ramirez-Ruiz 2004; Guetta et al. 2004; 
Pian et al. 2006) than luminous GRBs, which have been observed in large numbers out to 
higher redshifts. 
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Fig. 1.— Lightcurves of GRB 980425 detected with the brightest BATSE LAD 
(20-2000 keV; top) and with the BeppoSAX WFC (2-27 keV; bottom) plotted with 1-s res- 
olution. The four time intervals used in the current analysis are labeled with A, B, C, and 
D. The dotted lines indicate background levels. 
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Fig. 2. — Spectral parameter evolution of GRB 980425; the data points correspond to the 
intervals indicated in Figure 1. The uncertainties are la. 
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Fig. 3.— The HETE-2/FREGATE lightcurve of GRB 030329 generated using the publicly 
available Burst Lightcurve data integrated in 64-ms resolution time bins. 
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Fig. 4.— The INTEGRAL IBIS/ISGRI lightcurve of GRB 031203 in 1-s resolution. The four 
time intervals used in the current analysis are indicated with the dashed lines. The dotted 
lines show the background levels. 
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Fig. 5. — Single power law spectral index evolution of GRB 031203; the data points corre- 
spond to the intervals indicated in Figure 4. 
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Fig. 6.— The lightcurve of GRB 060218 seen with Swift BAT (15-350 keV; top) and XRT 
(0.2— lOkeV; bottom); the data are binned with 1.6 and 1.0s resolution, respectively. The 
14 time intervals used in our time-resolved analysis are indicated with dotted lines. 
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Fig. 7. — The spectral parameter evolution of GRB 060218. The last two points in the photon 
index plot (middle panel, marked with squares) were estimated by fitting a power law with 
exponential cutoff with the exponential cutoff energy fixed to the last well-determined value 
(4.7keV). All the uncertainties are la. 
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Fig. 8. — The evolution of the 7-ray isotropic-equivalent luminosity (2— 500 keV, source- 
frame energy) for all four events, in the source-frame time. For GRB 060218, the luminosity 
of the thermal component is plotted separately. 
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Fig. 9. — The unabsorbed best-fit duration-integrated spectra (solid lines) of the four events, 
overplotted with the deconvolved data (gray crosses) in the source-frame energy. The data 
are binned for display purposes. The analysis tool for the HETE-2 data of GRB 030329 was 
not publicaly available, and we only show here the spectral model presented in Vanderspek 
et al. (2004). See also text in § 2.5.1. 
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Fig. 10.— The £ pcak of the duration-integrated spectra of GRBs 980425, 030329, and 
GRB 060218; for GRB 031203, we plot a la lower limit. As a comparison, the £ pca k dis- 
tributions of 251 bright BATSE GRBs (Kaneko et al. 2006) and 37 HETE- 2 GRB/XRFs 
(Sakamoto et al. 2005) are plotted. Only well-constrained £ pea k values are included in the 
distributions. 
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Fig. II.— Locations of GRBs 980425, 031203, 030329, and 060218 in the £ pC ak-£ 7 ,iso plane. 
The dashed line indicates the correlation found by Amati et al. (2002). The events presented 
in Ghirlanda, Ghisellini, & Lazzati (2004) are also shown here. 
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Fig. 12. — The hardness ratio evolution of all four events in the source-frame time. The 
hardness ratio is determined in the source-frame energy for each event. The values for 
GRB 030329 were estimated from the spectral parameters presented in Vanderspek et al. 
(2004). 
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Fig. 13.— The evolution of the £ peak in GRBs 980425, 030329, and 060218, in the source- 
frame time. The dotted lines show the best-fit power law decay, £ pea k oc t^, with <fi = 
-1.46 ± 0.12, -1.17 ± 0.08, and -1.40 ± 0.06, respectively. For GRBs 980425 and 060218, 
the first points are excluded in the power-law fits. 
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Fig. 14. — Evolution of the X-ray (0.3— lOkeV, source-frame energy) prompt and afterglow 
luminosity (isotropic equivalent) in the source-frame time. 
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Fig. 15. — Cumulative total emitted energy (isotropic equivalent) of the four SN-GRBs 
(black lines) in 0.3— lOkeV (source-frame energy) as a function of source-frame time. As a 
comparison, the cumulative energy of the 10 Swift GRBs with known redshifts published in 
Nousek et al. (2006) are also shown here in color. 



-42- 



1Q 54 


III 


. □ □ b . 




□ ° 




^ 10 52 

CD 


□ 

• 

GRB 030329 


D n D D GRB 030329 
• 


o 

& 10 50 


□ 

GRB 031203 
• 

• GRB 060218 


- □ 

GRB 031203 
• 

GRB 060218 
• 


1Q 48 


GRB 980425 
• 

iii 


GRB 980425 
• 



10 47 10 48 10 49 10 50 10 51 10 52 0.1 1.0 10.0 100.0 1000.0 
Ex, iso (erg) T90,x (days) 



Fig. 16. — Comparison of the four SN-GRBs (full circles) with the 10 Swift GRBs presented 
by Nousek et al. (2006) (squares), in the E lt i so -Ex,i so plane (a) and in the -E^iso-Tgcx plane 
(b). The -E 7i i SO values used here are determined between 20— 2000 keV in the source frame. 
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Fig. 17. — Radio lightcurves for all four events at 4.8/4.9 GHz. The errors associated with 
the data are very small for GRB 980425 and GRB 030329. 
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Fig. 18. — Spline fits to the broadband radio lightcurves for the four events. 
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Fig. 19. — Cumulative isotropic-equivalent total emitted energy in 5—7 GHz (source-frame 
energy) as a function of source-frame time. 
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Fig. 20. — Isotropic-equivalent luminosity of SN-GRB X-ray afterglows scaled to t = 10 hr 
(source frame; 2—10 keV) after the burst trigger as a function of their isotropic 7-ray energy 
release (20—2000 keV). The redshift for each event is also shown in color (adopted from 
Nousek et al. 2006). 
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Fig. 21. — Cumulative isotropic-equivalent total emitted energy in 3000— 24000 A (source- 
frame energy) as a function of source-frame time, for all SNe associated with the GRBs 
presented in this work. 
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Fig. 22. — Summary of the isotropic-equivalent total emitted energy of the prompt and 
afterglow emission of the four GRBs, along with the properties of their associated SNe. 
The energetics in 7-ray, X-ray, Radio, and Optical (for SNe) wavelengths correspond to 
1-10,000 keV, 0.3-10 keV, 5-7 GHz, 3000-24000 A in the source frame, respectively. 



- 49 - 



1 1 1 1 1 
980425 * ' 030329 h_, 




a 


031203 < 1 

060218 < 1 

970828 ™ 22B « ' 

Q71 91 A . 
QROfil^. 1 . 

980703j_i 

990123 | | 






990510 i-i 

99070" 1 

99071 2 i i 

991216 | | 


990506 *- 




0002 1 4 < 1 

000 8162 , 2T^ H 


000131 * 1 




n 1 092 1 < 1 

01 121 1 i—i 

020124 | | 

020405 | | 






020813 <-^ nQ 

02*121 1 i i 

03022&I 1 

030328 i_i 
O^OA?Q , , 






030723 < 1 

1,1,1 






1 48 1 50 1 o 52 




1Q 54 


GRB E y (erg) 







199^bw m 



2003dh+ 
2003lw i— i— i 
2006aj + 
1997dq + 
1997ef + 
2002ap i— i— i 
2003jd + 

2005bf i-i-h 



SN E k (erg) 
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Table 1: Summary of the spectral fit results of GRB 980425. The data were fitted with a 



Comptonization model; all 


uncertaint 


ies are la. 






Time 




£peak 


a 


X 2 /dof 


Interval 




(keV) 






A 


17 A ±3.6 


175 ±13 


-0.12 ±0.22 


134.4/141 


B 


10.7 ±1.1 


133 ± 8 


-1.16 ±0.09 


149.9/141 



C 4.2 ±0.9 34 ± 3 -1.54 ±0.09 158.8/141 
D 2.0 T_ H 14 t I -1.51 ±0.36 121.8/141 



in units of 10 3 ph s 1 cm 2 keV 1 . 



Table 2: Summary of the spectral fit results of GRB 031203 using a single power-law model; 
all uncertainties are la. 



Time 


Amplitude 11 


Photon 


X 2 /dof 


Interval 




Index 




- 2 s 


6.56 ±2.46 


-1.60 ±0.09 


11.2/7 


2 - 4 s 


10.25 ±4.26 


-1.75 ±0.11 


4.9/7 


4 - 8 s 


4.81 ±1.66 


-1.67 ±0.09 


6.2/7 


8 - 50 s 


0.72 ±0.56 


-1.63 ±0.20 


5.0/7 



a in units of ph s 1 cm 2 keV 1 , at 1 keV. 
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Table 3: Summary of the joint BAT-XRT time-resolved spectral fits for GRB 060218. The 
spectra were fitted to a power law with high-energy cutoff with a blackbody unless otherwise 
noted. Note that we present £ pca k instead of a cutoff energy. The uncertainties are la. 



Time 


Nh 


kT 


^pcak 


Photon 


X 2 /dof Flux (0.5-150 keV) a 
(10~ 9 ergs -1 cm- 2 ) 


Interval 


(10 22 cm- 2 ) 


(keV) 


(keV) 


Index 



-8 ■ 


- 140 s b 










24.9 


±6.0 


-0 


.87 


±0. 


.75 


52 


.4/57 


6.38 ±4.00 


140 


- 303 s 


0. 


.47 


±0.08 


0.22 ±0.05 


36.1 


±7.2 


-1 


.39 


±0. 


.08 


238 


.3/292 


8.95 


+0.43 
-1.40 


304 


- 364 s 


0. 


.66 


±0.14 


0.17 ±0.03 


20.8 


±4.1 


-1 


.49 


±0. 


.11 


165 


.7/190 


12.84 


+0.49 
-2.87 


406 


- 496 s 


0. 


.62 


±0.11 


0.17 ±0.02 


31.5 


±5.0 


-1 


.39 


±0. 


.08 


236 


.4/302 


14.93 


+0.47 
-1.85 


496 


- 616 s 


0. 


.58 


±0.02 


0.16 ±0.02 


22.3 


±2.9 


-1 


.33 


±0. 


.06 


323 


.6/381 


13.99 


+0.47 
-1.12 


616 


- 736 s 


0. 


.51 


±0.07 


0.20 ±0.02 


18.9 


±2.6 


1 


.21 


±0. 


.07 


395 


.6/401 


12.59 


+0.39 
-1.02 


736 


- 856 s 


0. 


.64 


±0.08 


0.16 ±0.02 


15.2 


±2.4 


1 


.41 


±0. 


.07 


397 


.1/412 


13.40 


+0.37 
-1.04 


856 


- 976 s 


0. 


.66 


±0.08 


0.15 ±0.01 


12.0 


±2.0 


-1 


.45 


±0. 


.07 


390 


.2/420 


10.03 


+0.39 
-1.37 


976 


- 1256 s 


0. 


.61 


±0.06 


0.15 ±0.01 


5.7 


±0.4 


-1 


.39 


±0. 


.09 


611 


.7/567 


11.17 


+0.23 
-0.93 


1256 


- 1556 s 


0. 


.65 


±0.06 


0.14 ±0.01 


3.6 


±0.2 


-1 


.30 


±0. 


.14 


571 


.2/539 


10.26 


+0.30 
-0.97 


1557 


- 1857 s 


0. 


.07 


±0.06 


0.13 ±0.01 


2.4 


±0.4 


-1 


.47 


±0. 


.03 


483 


.5/481 


8.61 


+0.39 
-0.95 


1857 


- 2157 s 





.60 


Fixed 


0.14 ±0.003 


1.8 


±0.7 


-1 


.62 


±0. 


.15 


427 


.8/433 


6.41 


+0.49 
-0.62 


2157 


- 2457 s c 





.00 


Fixed 


0.14 ±0.004 






-2 


.45 


±0. 


.03 


394 


.8/388 


6.34 


+0.42 
-0.27 


2457 


- 2734 s c 





.60 


Fixed 


0.13 ±0.004 






-2 


.54 


±0. 


.04 


327 


.1/342 


6.15 


+0.37 
-0.52 



a Unabsorbed flux. Uncertainties are associated with the absorbed flux estimated from the 

fitted parameters. 

b Only BAT Event data were used. 

c Fitted by a power law. 
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Table 4: Summary of broadband properties of the prompt and afterglow emission of the 
four GRBs, along with the properties of their associated SNe. The energetics in 7-ray, 
X-ray, Radio, and Optical (for SNe) wavelengths correspond to 1 — 10, 000 keV, 0.3— lOkeV, 
5-7 GHz, 3000-24000 A in the source frame, respectively, unless noted. 





980425 


030329 


031203 


060218 




(1998bw) 


(2003dh) 


(20031w) 


(2006aj) 


z 


0.0085 


0.1685 


0.105 


0.0335 


T 90, 7 (s) 


34.9 ±3.8 b 


22.9 C 


37.0 ±1.3 


2100 ±100 d 


Sx,2-30 kcV a 


1.99E-6 


6.71E-5 


8.46E-7 


1.09E-5 


5*7,30-400 keV a 


3.41E-6 


1.20E-4 


1.74E-6 


3.09E-6 


Sx/ Sy 


0.58 


0.56 


0.49 


3.54 



£ 7jiso (erg) 9.29 ±0.35 xlO 47 1.33 xlO 52 1.67 lg;^ xlO 50 4.33 +?;^ xlO 49 
£ pcak (keV) 122 ±17 70 ±2 C > 71 4.7 ±1.2 



T 9 o,x (day) 


640.1 


2.61 


79.8 


9.31 


£x,iso (erg) 


1.67 xlO 48 


7.09 xlO 50 


8.27 xlO 48 


6.15 xlO 47 


T 90 ,r (day) 


68.6 


72.4 


81.2 


67.2 


Er; 1so (erg) 


4.21 xlO 44 


5.64 xlO 46 


1.41 xlO 45 


1.09 xlO 43 


T 90 ,sn (day) 


53.9 


31.4 


53.6 


30.7 


^SN.iso (erg) 


2.31 xlO 49 


1.81 xlO 49 


3.15 xlO 49 


9.24 xlO 48 



a in units of erg cm~ 2 . 

b Reference: BATSE current catalog (http://gammaray.nsstc.nasa.gov/batse/grb/catalog/current/). 
c Reference: Vanderspek et al. (2004). 
d Reference: Campana et al. (2006). 
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Table 5: Log of the radio observations for GRB 060218. The flux calibration has been per- 
formed using the source 3C286 as standard. The observations on February 21 and February 
28 were alternated between 1.4 and 4.9 GHz, and 2.3 and 4.9 GHz, respectively, observing in 
40 minute blocks at each frequency. 



Date (2006) 


AT 


Integration Time 


Frequency 


Flux 




(days since trigger) 


(Hours) 


(GHz) 


(A*Jy) 


Feb 21.452 - 21.951 


3.30 - 3.80 


5.9 


1.4 


89 ± 108 


Feb 28.463 - 28.908 


10.31 - 10.76 


5.2 


2.3 


-35 ± 70 


Feb 21.481 - 21.927 


3.33 - 3.78 


5.2 


4.9 


257 ± 36 


Feb 24.688 - 24.943 


6.54 - 6.79 


5.9 


4.9 


125 ± 38 


Feb 28.433 - 28.932 


10.28 - 10.78 


5.8 


4.9 


106 ± 42 


Mar 10.642 - 10.905 


20.49 - 20.76 


6.1 


4.9 


6 ± 34 


Apr 1.345 - 1.844 


42.20 - 42.70 


12.0 


4.9 


8 ± 24 
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Table 6: Lower limits on the Lorentz factor of the material responsible for the prompt 
emission (following the formalism of Lithwick & Sari 2001). The high-energy photon index 
(f3) and £ pea k values are for the peak spectra of duration 5T, where the parameter limits are 
la. 



GRB 




5T 




P 




£ b 

'--max 




Limit 


Limit 






(s) 




(keV) 


(keV) 


A 


B 


980425 


0.0085 


12.6 


0.00314 


>3.5 


143 


300 


330 


1.6 


2.4 


030329 


0.1685 


5.27 


0.70 


2.44 


52.5 


400 


4.9 x 10 7 


13 


27 


031203 


0.105 


8.16 


<0.002 


2 
3 


>148 


200 


1.0 x 10 6 
1.5 x 10 5 


8.9 
3.7 


16 
7.5 


060218 


0.0335 


1981 


0.020 


2.75 


8.8 


100 


0.36 







a The flux normalization at £ pca k, F p = f£~f ak , in units of ph s -1 cm -2 keV 
b The maximum photon energy with significant detection. 
c The optical depth for a photon of energy m e c 2 for 7/3 « 1. 



Note. — For f < 1 and £ max < m e c 2 (as is the case for GRB 060218) there is no lower limit on the 
Lorentz factor 7. Limit B derivation assumes that the high-energy power law continues significantly above 
the observed range for the cases presented in the Table. For GRB 031203 there is only a lower limit on £ pea k, 
and (3 is not measured, so we conservatively take £ pea k to be equal to its lower limit and show results for 
two representative values of (3. 



Table 7: The source-frame time, tx, of the peak in the £Lx,iso history and the corresponding 
^x-^x,iso(^x) values for the four events. 



GRB 


z 


tx 
(s) 


^X-^X,iso(^x) 

(erg) 


txLx,iso(tx) 
/-^7,iso 


980425 


0.0085 


2.4 x 10 7 


6.3 x 10 47 


0.68 


030329 


0.1685 


1.7 x 10 4 


5.3 x 10 50 


0.040 


031203 


0.105 


3.5 x 10 5 


2.1 x 10 48 


0.013 


060218 


0.0335 


3.4 x 10 4 


2.9 x 10 47 


0.0068 
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Table 8: Summary of the estimated minimal combined energy, E min , for all four events, with 
v R = 4.86 GHz. 



GRB 


z 


(days) 


(mJy) 


E ■ 

- L 'min 

(erg) 


(cm -3 ) 


980425 


0.0085 


200 


3 


1.5 x 10 48 


0.0022 


030329 


0.1685 


100 


3 


1.8 x 10 49 


0.32 


031203 


0.105 


100 


0.35 


3.1 x 10 48 


0.048 


060218 


0.0335 


7 


0.16 


2.1 x 10 46 


0.78 



